“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1996-03 


Evaluation of the Boeing PAN AIR technologies 
code (A502i) through prediction of separation 
forces on the GBU-24 


LeTourneau, Matthew A. 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/8328 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


| (8 D U DLEY research materials and institutional publications created by the NPS community. 
. | Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
ith 
KNOX appointed — and published — scholarly author. 


INN LIBRARY Dudley Knox Library / Naval Postgraduate School 


411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 





THESIS 


EVALUATION OF THE BOEING PAN AIR 
TECHNOLOGIES CODE (A502I1) THROUGH 
PREDICTION OF SEPARATION FORCES 
ON THE GBU-24 
by 


Matthew A. LeTourneau 


March, 1996 


Thesis Advisor: Max F. Platzer 





Approved for public release; distribution is unlimited. 


Thesis 
L55513 


OUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-51 01 


—— ar 
Form Approved OMB No. 0704-0188 " 





REPORT DOCUMENTATION PAGE 


Public reporung burden for this collecuon of information ts estimated to average | hour per response, including the time for reviewing instruction. searching existing 
| data sources, gathering and maintaining the data needed. and completing and reviewing the collection of information. Send comments regarding this burden estimate 
or any other aspect of Unis collection of information, including suggestioas for reducing this burden, tu Washington (lcadquarters Services, Directorate for lifurmation 
Operations and Reports. 1215 Jefferson Davis Highway. Suite 1204. Arlington, VA 22202-4302. and to the Office of Management and Budget. Paperwork 


Reduction Project (0704-0188) Washington DC 20503. | 


AGENCY USEONLY (Leave blank) aa REPORT DATE . REPORT TY PEANDDATES COVERED 
March 1996 Master's Thesis 


4. TITLEANDSUBTITLE EVALUATION OF THE BOEING PAN AIR 5S. FUNDINGNUMBERS 
TECHNOLOGIES CODE (A502I) THROUGH PREDICTION OF 
SEPARATION FORCES ON THE GBU-24 


6. AUTHOR(S) LeTourneau, Matthew, A 


7. PERFORMINGORGANIZATIONNAME(S) AND ADDRESS(ES) 8. PERFORMINGORGANIZATION 
Naval Postgraduate School REPORT NUMBER 
Monterey CA 93943-5000 


9. SPONSORING MONITORING AGENCY NAME(S) ANDADDRESS(ES) 10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the 
official policy or position of the Department of Defense or the U.S. Government. 







12a. DISTRIBUTION: AV AILABILITY STATEMENT 


12b. DISTRIBUTIONCODE 





> 
4 
os 
t 


Approved for public release; distribution is unlimited. 


13. ABSTRACT (maximum 200 words) 
The Boeing PAN AIR Technologies code (A5021) is investigated to explore its suitability for 


determination of separation forces on ordnance. To this end, A5021 is first assessed by applying it to 
three problems for which other solutions and experimental data are available, 1.e. steady flow past a 


rectangular, parabolic arc wing and a delta wing at both subsonic and supersonic conditions. Good 
avreement is found in all vases. A502i is then applied to the GBU-24’s being in two configurations for 
a subsonic case and a supersonic case. Good agreement is found with data obtained from wind tunnel 





experiments for low angles of attack. 


14. SUBJECT TERMS A502i, Panel Methods, Potential FLow, Geometry Modelling, 
GBU-24 


15. NUMBEROF 
PAGES 116 
16. PRICECODE 


20. LIMITATIONOF 





















17. SECURITY CLASSIFICA- | 18. SECURITY CLASSIFI- 19. SECURITY CLASSIFICA- 
TION OF REPORT CATION OF THIS PAGE TION OF ABSTRACT ABSTRACT 


Unclassified Unclassified Unclassified UL 


NSN 7540-01 -280-5500 Standard Form 298 (Rev 2-89) 
Prescribed by ANSI Std. 239-18 298-102 












Approved for public release; distribution is unlimited. 


EVALUATION OF THE BOEING PAN AIR TECHNOLOGIES CODE (A5021) 
THROUGH PREDICTION OF SEPARATION FORCES ON THE GBU-24 


Matthew A. LeTourneau 


Lieutenant, United States Navy 
B.S., The George Washington University, 1988 


Submitted in partial fulfillment 
of the requirements for the degree of 


MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 
from the 


NAVAL POSTGRADUATE SCHOOL 
March 1996 





DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOO). 
MONTEREY CA 93943-5101 


ABSTRACT 


The Boeing PAN AIR Technologies code (A502i) is investigated to explore its 
suitability for determination of separation forces on ordnance. To this end, A5021 is first 
assessed by applying it to three problems for which other solutions and experimental data 
are available, i.e. steady flow past a rectangular, parabolic arc wing and a delta wing at 
both subsonic and supersonic conditions. Good agreement is found in all cases. A5021 is 
then applied to the GBU-24’s being in two configurations for a subsonic case and a 
supersonic case. Good agreement is found with data obtained from wind tunnel 


experiments for low angles of attack. 
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I. INTRODUCTION 


In the past, ballistic trajectory determination for manual or computer predicted 
ordnance delivery from an aircraft was determined through measurement of separation 
forces on the piece of ordnance via wind tunnel or captive carry measurements. The advent 
of panel method codes using linearized potential theory, such as A5021, or its full potential 
version TranAir, allow for a cheaper and safer method of predicting separation forces. 
Furthermore, A502i allows for any arbitrary configuration to be modelled within the 
limitations of the number of panels and networks allowed and excluding transonic flow. 

The purpose of this work was to determine the separation forces on a GBU-24 
carried by an F-14 on stations 3 or 6 or both. It was also the purpose of this work to 
provide an analysis of the code itself to see if it is a viable tool for the study of flow 
characteristics over arbitrary wing configurations for use in the Naval Postgraduate 
School’s (NPS) Department of Aeronautics and Astronautics. The majority of the work 
was conducted on the NPS computer systems. The Department of Aeronautics and 
Astronautics Silicon Graphics Incorporated (SGI) workstations were utilized for most of 
the input files as well as the execution of the code. Due to the amount of disk space 
required, storage of the output files took place on the NPS Computer Center’s Y-MP EL98 
Cray computer. The bulk of the GBU-24 data was calculated using the SGI workstations 
at the Naval Air Warfare Center in Warminster. 

The scope of this analysis was to understand the capabilities of the A5021 code. 
The approach was to validate A5021 against existing data and linear theory. The code was 
run for three different geometries under assorted Mach and AOA conditions. Comparisons 


were made for each of the geometries. 








Il. OVERVIEW OF THE A502i CODE 


The A5021 code is used to computationally analyze inviscid subsonic or supersonic 
flows about arbitrary configurations. The code differs from other panel methods in that it 
is a higher order panel method; that is, the singularity strengths are not constant on each 
panel. A502i solves the linearized potential flow boundary-value problem at subsonic and 
supersonic Mach numbers. 

The aerodynamic solution provides surface flow properties (flow directions, 
pressures, Mach number), configuration forces and moments, sectional forces and 
moments, and pressures. Additionally, A502 calculates flow properties in the flow-field 
points and flow-field streamlines. Results are limited to subsonic and supersonic cases 
(transonic cases excluded) with attached flow. Results are not usually applicable to cases 


where viscous effects and separation are dominant. 
A. THEORY 
The basic equations describing the flow of a viscous, compressible, heat- 


conducting fluid are the Navier-Stokes equations. These are: 


(a) The continuity equation, 
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X =(X,,X.,X,), and where the conventional index notation is used instead of X=(x,y,z). In 
addition, tis time, p(X,t) is the density and V(X,t) is the velocity vector, with components 


V=(V,,V.,V;). 


(b) The conservation of momentum equation, 
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Where t,, is the deviatoric portion of the stress tensor which vanishes for a frictionless 
fluid, f(X,t) is an external body force per unit mass exerted on the fluid, and p(X,t) is the 


pressure. 


(c) The conservation of energy equation 
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where e(X,t) is the internal energy of the fluid, k 1s the coefficient of heat conductivity for 


the fluid, and T(X,t) is the temperature. 
(d) The equation of state 
f(p, p, P)=0 (2.4) 
where the function f depends on the type of fluid. The derivations of these equations can 
be found in Refs. 1,2 and 3. 
The Navier-Stokes equations can be simplified by the neglect of viscosity, which is 
equivalent to setting the deviatoric stress tensor to zero. Combining the momentum and 


continuity equations yields 
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where the convective derivative operator is defined as 
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Equation (2.5) is referred to as Euler’s equation. The continuity and energy equations can 


Equation (2.5) is referred to as Euler's equation. The continuity and energy equations can 


be reduced to 
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and the rate of increase of heat per unit mass 1s given by 
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Equations 2.5, 2.6 and 2.7 can be reduced to a single equation when four further 
assumptions are made. First, assume isentropic flow, thus 
q=0 (2.8) 
Second, assume irrotationality 


VxV=0 (2.9) 


which allows for the introduction of the potential function [Refs 2,3] 


V@=V (2.10) 


Third, assume the existence of a freestream potential ®_, whose gradient is the uniform 


velocity V_ attained at points sufficiently distant from the disturbance being analyzed, and 


thus write 
o=0-® a 


and 


V=(u,v,w)=Vd=VO +VO=V_+V0 By) 
The quantities @ and V are called the perturbation potential and velocity [Ref 3]. Fourth, 


assume that 
WP <a, (2. 13) 


everywhere, where a_ is the freestream speed of sound. 


Based on these four assumptions, the unsteady potential equation is obtained 


[Refs 1.3]: 
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Assuming the flow conditions do not change with time yields the steady non-linear 

potential equation. 
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(226s) 
When M_=0, equation (2.15) reduces to Laplace’s equation, 
V°o=0 (2.16) 
For the case of M_ # 0, the following is supposed, 
M?|¥|<1-M- (2.17) 
M:|¥\<1 (2.18) 


which are small perturbation assumptions [Refs. 1,2]. With these assumptions the steady 


non-linear potential equation reduces to the Prandtl-Glauert equation [Ref 1]: 


(1-M-)0,,+ >, +, =0 (2.19) 


Through a coordinate transformation [Refs. 1,2,3], the Prandtl-Glauert equation can be re- 
written as: 

SQ. t 955+ 0,=0 (2520) 
where when s=1, it is the subsonic case and Laplace’s equation applies and when s=-1, it 
1s the supersonic case and the wave equation applies. Applying Green’s third identity [Ref. 


l yields the following integral equation, 
An S R a R ( : ) 


where O represents the source strength and uw represents the doublet strength. When 


supplemented with boundary conditions, it 1s equation (2.21) that A502: solves. 
A502i solves equation (2.21) through a discretization process. The general idea of 
the process falls into two parts. The first is developing finite dimensional approximate 


representation formulas for the singularity functions, which creates a system of equations 


with unknown coefficients, A,. The second part involves solving the set of equations for 
allA.. This allows for completely determining the source and doublet functions. Then, by 


virtue of equation (2.21), the potential function o(P) is determined for all points P, solving 


the problem. 

The features of A502i which distinguish it from predecessors are three-fold. The 
first is a feature known as “continuous geometry”, the second ts linear source and quadratic 
doublet variation, the third is continuity of doublet strength. 

Most panel methods approximate the configuration geometry with panels whose 
planform is a quadrilateral. Thus, if the panels themselves are planar, only a small class of 
configurations (such as cylinders and flat wings) can be described without gaps being left 


between panels. These gaps are generally small, except for highly twisted surfaces. The 
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gaps cause little numerical error in subsonic flow, but in supersonic flow, the cumulative 
effect of the gaps is serious [Ref. 1]. The problem is not associated with leakage of flow 
through the gaps, but with the doublet strength jumping abruptly from a non-zero value to 
zero at a panel edge which does not exactly meet the adjacent edge. In A502i, gaps are 
closed by means of panels which are comprised of several planar regions. 

The feature of linear source and quadratic doublet variation is what makes A502i a 
higher order panel method. The basis function corresponding to a source parameter Is 
locally linear, while the basis function corresponding to a doublet parameter is locally 
quadratic. This is what allows for A5021 to find supersonic solutions. Numerical solution 
of the wave equation Is far more sensitive to the numerical idiosyncracies of a panel method 
than is the solution of Laplace’s equation for subsonic flow. Experimental evidence [Ref. 
1] indicates that exact surface analysis is not feasible in supersonic flow without doublet 
continuity, thus the potential for numerical error is greatly reduced by requiring the doublet 


singularity strength to be continuous across panels. 
B. GENERAL A5021 USAGE 


The use of the A502i code consists of generating an input file, which can be 
arbitrarily named, and which contains the information defining the geometry of the 
configuration, flow field points of interest, the flow conditions and wakes. The process of 
building a geometry is difficult in that A502i is particular about its input format. Simple 
configurations, such as a rectangular, planform wing can be modelled manually, but more 
complex structures require a pre-processing program, such as MACGS, where a geometry 
can be graphically built. MACGS will output a data file in a format that, with minor 
modifications, via another pre-processing program that can move the data from three 
columns to six columns, will be readily usable by A502i. Currently, the school does not 


have a copy of MACGS, but it can be acquired through McDonnell-Douglas. To complete 


this thesis, MACGS was used on the SGI workstations at NAWC Warminster. Wakes 
also must be constructed in the same manner as the structure to be analyzed. More detailed 
instructions on the specifics of wakes and surface geometries can be found in Ref. [4]. 
Appendix A is a portion of an output file, but lines 1 thru 1120 are an exact duplicate of the 


input file. 
1. Running A502i with an Existing Executable 


Assuming an A5021 executable file (e.g., A502) has already been placed in a user 
executable directory (e.g., /usr/local/bin), the only other necessary items needed to produce 
a set of A5021 output files is the input file and a large amount of storage space. Anything 
modelled with more than one thousand panels total will use more than one hundred mb of 
disk space. If the maximum number of panels (20,000) is used, the disk space required 
will be on the order of 2 gb. 

To run A5021, enter after the UNIX prompt: 

A502 <input file > output file 
Prior to running the code, it is highly recommended that a Cray account be opened and 
linked to the department’s SGI workstations. This is done by assigning the same user i.d. 
number to the Cray account as is assigned to the account with the department. User i.d. 
numbers can be changed by the computer center at the user’s request. This is required due 
to the limited disk space available to individual accounts in the department. Once an 
account is opened, log on to a department workstation, change directories to an existing 
Cray directory, for example (after the UNIX prompt): 
cd /jedi/d1/maletour 

Transfer the input file to the Cray directory and execute the code. The screen will display 
what portion of the code it is performing and how long it took to perform each portion in 


CPU time. The code outputs numerous files in addition to the arbitrarily named output file. 


The two output files of interest, in the vast majority of cases, are the arbitrarily named 
output file and the ft13 file. In order to run another solution all output files must be deleted 
or renamed prior to re-executing the code. Relevant results should not be kept on the Cray 
account as files on disks dl, d2 and ul are considered temporary storage and subject to 


erasure after a period of time. 
2. Creating an A502i Input File 


The input file, which can be arbitrarily named, consists of two portions, the largest 
being the geometry data. Appendix A is a complete recreation of the input file for the 
GBU-24 with canards. The file begins with line 1, $STITLE, and ends with line 1120, 
SEND. The line numbers are for reference only and are not part of the actual input file. 
The first portion consists of creating the initial conditions, i.e., the free-stream Mach 
number and angle of attack, the type of analysis to be performed, i.e., solution or 
datacheck, what types of output that are to be included in the arbitrarily named output file, 
and reference points to be used in calculating forces and moments. The geometry data 
consists of the points that bound each panel, that in turn belong to a specific group of 
panels that make up a network. The overall structure being modelled consists of a series of 
networks. A50O2i can run up to 150 networks and or 20,000 panels with a limit of 8,000 
panels per network. Referencing Appendix A, line 28 represents the first network of the 
model, a canard. Line 29 represents the number of networks that will be classified under 
this S$POINTS statement. Line 30 indicates what kind of surface the network will be, a 
three-dimensional surface with flow properties to be calculated, a wake and a base are 
several examples. Line 31 is the number of y points and the number of x points 
respectively that make up the grid of that network. Line 32 is where the panel points start. 
Reference 4 contains detailed instructions on the options and meaning of each of the non- 


geometry inputs, including some capabilities not shown in Appendix A. 
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Two types of solutions can be run, a datacheck and a full solution. Reference 4 
explains how to enter either one into the input file. The datacheck only analyzes the 
geometry. This can be accomplished in a matter of seconds for a simple geometry as it is 
only running the first several portions of the code. The full solution can take a couple of 
hours for a geometry of the size of 4,000 panels. The datacheck should be run once the 
seometry has been modelled. It will check for any panel edges that do not abut properly, 
and when column 4 of line 20 in Appendix A is a 1, the datacheck will list the unit normal 
vectors, which must be facing outward. The datacheck will also see if the wakes are 
attached properly. A5021 1s capable of giving warnings both on-screen and in the 
arbitrarily named output file when an edge or a wake 1s not modelled properly, but it only 
lists the unit normal vectors. The directions of the vectors must be manually checked by 
the user. The full solution performs the datacheck first, so the data is repeated in the 
arbitrarily named output file. Appendix B is a portion of the output file that contains the 
summary of facing surfaces. Each panel edge is looked at to see what it abuts against. 
Sections such as wingtips, leading edges of a flat plate or any surface that does not need a 
wake attached, but is unabutted to any other panel on that edge will draw probable error 
messages or warnings from the code. The user must ensure that the edge is not supposed 
to abut against anything or need a wake attached. If that is the case, the warnings may be 
ignored. Appendix D is the first page of the portion of the output file that lists the unit 
normal vectors. The three columns under zc are the x-y-z coordinates of the given panel’s 
center. The three columns under znc are the x-y-z coordinates of the unit normal vector. In 
most cases, when the y coordinates are of the same sign, then the unit normal vector is 


pointing outward. 
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C. GEQMETRY MODELLING 


Five geometries needed to be modelled, each of increasing complexity. Modelling 
proved to be the most difficult task. in that A502i is a FORTRAN code and is very format 
sensitive, but the sheer number of points that need to be generated can take a lot of time and 
the order those points are listed in the input file is what determines whether or not the shape 
is correctly modelled. Of the five geometries modelled, none were done completely 
manually. A spreadsheet was used for generating the parabolic arc airfoil and the deltawing 
since those structures can be constructed out of one network, excluding wingtips and 
wakes, and the surface can be defined by a mathematical function. The bombs and the F- 
14 required the use of MACGS to be properly modelled. MACGS is indifferent as to the 
order that geometries are built, and often doesn’t require many coordinate inputs if building 
a model on top of an existing IGS file. The output file from MACGS is automatically 
formatted and the points placed in the appropriate order for A502i to understand. 
Although, the order may be reversed where the unit outward normal vector is concerned. 
MACGS has the ability to output files in several different panel method code input formats, 
including PMARC. Reference 4 gives detailed instructions on how to properly order 
points to build a group of networks that will model a geometry. A5021 uses a right-handed 
coordinate system that is similar to an aircraft body axes. When put in terms of a wing, the 
X axis is positive from leading edge to trailing edge. The z axis is positive up and the y axis 
is positive out the night wing 


1. The Parabolic Arc Airfoil 


The parabolic arc airfoil is the simplest of all the geometries. The airfoil has a chord 
of five and a span of ten. The maximum thickness is .15. The model consists of 
approximately 600 panels, including the wake and wingtips. A spreadsheet was used to 


develop the geometry portion of the input file. Line 32 of Appendix A demonstrates the 


RZ 


format that the spreadsheet used. Rows consist of two points, with coordinates x1, yl, z1, 
x2, y2, z2 using a format of 6F10.0. The chord was divided into 25 points (x coordinate) 
from trailing edge to leading edge and then another 25 points from leading edge to trailing 
edge (bottom half). The span was divided into 12 points (y coordinate) from left to right. 
Due to the symmetry of a rectangular planform, the y coordinate was constant along the 50 
x coordinates that constituted a chordwise cross-section. To attain a maximum thickness of 


.15 the formula, 





z=.34(—- =) (2.22) 


was utilized to generate the values of the z coordinates. The wingtips simply connect the x 
coordinate on the top side with it’s symmetrical counterpart on the bottom side. Due to a 
trailing edge composed of a straight line, the wake is modelled by a single panel that spans 
the trailing edge and has a length aft of 100. Figure 2.1 shows the panel distribution across 
the top surface of the parabolic arc airfoil, where the thickness is represented by the color 


scheme. A panel and a point are numbered to show how they were entered into the input 


file. 


2. The Deltawing 


The deltawing represented a step up in complexity over the parabolic arc airfoil. 
The chordwise cross-section is parabolic, while the spanwise cross-section is linear. The 
procedure for building the geometry on a spreadsheet was the same as that for the parabolic 
arc airfoil, only the chord length is not constant along the span. For simplicity in design, 
the number of panels per column of panels is constant on the deltawing, as on the parabolic 
arc airfoil. This means an increasing panel density in the direction of the wing tip. The 
wake is modelled the same as the airfoil. The right wingtip ended in a point, so no extra 


panelling was needed to close any gaps. The symmetry of the deltawing allowed for 
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Figure 2.1 Parabolic Arc Panel Distribution 
further simplification and reduction of the code’s run time by only modelling from the 
centerline to the right tip. A502i allows the user to stipulate whether there is symmetry in 
the x-z plane and or the x-y plane (see line 5 and 6 of Appendix A). This means that the 
gap between the top and bottom panels at the center line does not need to be bridged as in 
the parabolic arc airfoil (symmetry could have also been used in the airfoil’s case). The 
chord of the deltawing has a length of 90 and the semi-span has a length of 15. The 
maximum thickness occurs midway along the centerline and is .05. The model consists of 


880 panels. Figure 2.2 shows the panel distribution along the top surface of the deltawing. 





Thickness is represented by the color scheme. 
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Figure 2.2 Panel Distribution of Deltawing 
3. GBU-24 


Wind tunnel experiments were run on the GBU-24 without canards attached, so it 
was deemed relevant to build a model with and without canards as a comparison of the 
code’s performance. The model of the GBU-24 was too complex to build with a 
spreadsheet, so MACGS was used. The bomb was modelled at NAWC Warminster by 
superimposing a group of networks on top of an IGS file being displayed by MACGS. 
The complete configuration consists of approximately 1300 panels. Figures 2.3 and 2.4 


are displays of the GBU-24 with canards, with Figure 2.4 including the wakes. Figure 2.5 


IS 


is included to show how the GBU-24 model was assembled. Each different color 
represents a network. 

Several features of the geometry are relevant to point out. Two of them are 
modifications made to the geometry that differ from the actual dimensions of the bomb. 
Dr. Alex Cenko of NAWC Warminster has extensive experience with modelling stores in 
A5021. The modifications were made on his knowledge of how to get the most accurate 
results from the code when modelling stores. The first is to model the fins and canards as 
flat plates, i.e., no thickness, which A5021 allows you to do through a single numerical 
change in the input code for each network that represents a flat plate (see line 30 of 
Appendix A). The fins and canards are extremely thin when compared to the rest of the 
bomb, and to add a third dimension to the geometry complicates the construction of the fin 
or canard for several reasons. The leading and trailing edges must be sharp and the surface 
the fin or canard attaches to would have to be modified to abut properly with two edges 
instead of one. Experience has shown that the simpler version yields accurate predictions. 
A502i is an inviscid code, so it cannot take into account Separation effects on its own. The 
GBU-24 does not have a flat base. In reality, it is more bullet nosed in shape. However, 
at the speeds with which the bomb 1s being analyzed, separation does occur near the trailing 
edge of the bomb. Experience has shown that truncating the end into a flat base and 
designating it a separated flow region through an appropriate input (see line 702, column 1, 
Appendix A) yields better results than attempting to model the bomb to exact physical 
dimensions. The last feature to point out are the wakes, as seen in Figure 2.4. A502i has a 
limitation in that the wakes must be modelled by the user, and they have the same abutment 
requirements as physical surfaces. Regardless of angle of attack, the wakes remain 
stationary with respect to the body to which they are attached. At higher angles of attack, 
the wakes are no longer close to paralleling the free-stream velocity. Remodelling the 


wakes is nearly an impossible task. The fin and base wakes would not be too difficult 
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Figure 2.4 GBU-24 Geometry (Wakes Shown) 
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because they do not abut against any physical surfaces except for the surfaces they are 
trailing from. The canard wakes must abut properly against the bomb’s body all the way to 
the base. Modification of those wakes would entail modification of the entire body, or ina 
possible simplification, letting the wakes remain attached along the body until the base of 


the bomb and then shifting them relative to the free-stream. 
4. The F-14 


The F-14 geometry was modelled in the same fashion as the GBU-24. The 
geometry consists of approximately 1500 panels. While that may seem fairly coarse for 
such a complex structure, experience shows that it is all that is required to get accurate 
predictions. The primary area of interest is the underside of the fuselage forward and 
between the two nacelles. Higher panel density on the top half is not required. Figures 2.6 
and 2.7 display the F-14 geometry without and with wakes shown. Several omissions are 
made to the model as having a trivial effect on the analysis or no effect at all. Phoenix rails 
and bomb racks are not modelled along with the chin pod because they are deemed 
insignificant to achieve reasonably accurate predictions over small angles of attack. 
External tanks were not considered, but could be modelled much in the same way as the 
bomb and inserted into the input file to see what effects the drop tanks have on separation 
forces. The vertical tails and horizontal stabilizers were deemed irrelevant to the prediction 
of the separation forces and were left out. This reduces the number of panels and 


networks, which also reduces the amount of time it takes to run a solution. 
5. Combination Geometries 


The F-14 and GBU-24 were modelled separately, but were combined together as 
shown in Figure 2.8. The first step to accomplish this was using the FORTRAN code 
NAVSEP which, among many of its functions will translate coordinates to relocate items in 


the flow-field. Once accomplished, the GBU-24 file was pasted into the F-14 input file. 
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Figure 2.8 F-14 and GBU-24 on Station 3 Combined Geometries 





D. GRAPHICS VISUALIZATION 


One of the pre/post-processing codes that came with A5021 is called RAID. 
Currently an executable exists in the department’s computer system that can be accessed 
through typing after the UNIX prompt, 

raid 
Raid is a basic graphics program that can read A5021 and TRANAIR input files and display 
geometries and flow properties from solutions. After accessing raid, it will ask what type 
of input file it is being asked to visualize. It can handle five other modes of input besides 
A502i/TRANAIR [Ref. 10]. It will then ask for the name of the input file. A prompt will 
follow asking about object definition matrices which needs to be answered by, 
EACH 

The next prompt will ask if the panels are going to be shaded by a Cp value. Cp is used 
generically in that Cp can be displayed or Mach or any of the relevant 49 surface flow 
properties [Ref. 4]. This 1s only used after a solution has been run, data has been extracted 
and a colorscale for contour plots has been determined. If only the geometry is to be 
displayed, then hit carnage return for all the next questions until a pink window appears 
with a menu in the lower left corner of the window. An anomaly of the program is that if 
you want to display the wakes, then wake display must be deselected. To view the 
geometry select view on the menu. All selections with the mouse in RAID are made with 
the center mouse button. A new window appears with a menu bar at the bottom and left 
and the geometry in the center. From there, rotation, translation, scaling, axes, reflections 
and other manipulation of the geometry is possible. Figures 2.3 through 2.8 are examples 
of geometries displayed on RAID. 

When presentation of solutions (i.e., Cp or Mach contours) are desired, the use of 


another post-processing program is required to generate a colormap file. The program is 
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called crebar. An executable currently exists on the department's system. Type in, 
crebar 

after the UNIX prompt. The program will ask straightforward questions. Number values 
associated with colors available can be found in Ref. 10. The color file can be saved under 
any name, but must lie in the same directory as the input and solution files. The first line of 
the color file will list four numbers. The first number 1s the number of colors assigned to 
the colormap (248 maximum). Occasionally, the color bar displayed in RAID when using 
a colormap will disappear when certain menu items are selected. To prevent this, change 
the last three numbers to read 6, I, -1. Plotting outputs from RAID requires saving the file 
in a format, such as RGB, that a printer will recognize. It is possible to change the text 
color and background color, the default is black, to avoid excessive use of black ink in hard 


copies. 
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WI. DATA EXTRACTION 


When a solution has run to completion, there are two files of interest, the arbitrarily 
named output file and the ft13 file. The arbitrarily named output file contains results for 
everything that A5021 solves for. The ft13 file contains only the 49 surface flow properties 
on each panel. Appendix Dis the solution portion of an arbitrarily named output file for the 
first network. 

For purposes of displaying flow properties on RAID, it is necessary to utilize the 
ftl3 file. A post-processing code called RAIDCONV is used to extract the specific 
information. To access RAIDCONV, type 

raidconv 
after the UNIX prompt. An executable currently exists in the department’s system. The 
ftl3 file must be in the same directory as RAIDCONV 1s accessed in. RAIDCONV will 
prompt the user for which kind of panel method 1s being used (A502i is one of three 
choices). The next prompt will ask for the name of the ft13 file. The last prompt will ask 
for the flow property that 1s to be extracted. A file called ftl3.cp will be created. It can be 
renamed for purposes of multiple flow properties extraction. Abbreviations for the 49 flow 
properties can be found at the bottom of page 1 of Appendix D. The two primary flow 
properties are 
LMACHU for local Mach number 
CP2ND for second order pressure coefficient 

CP2ND is the default setting for RAIDCONV. Once the ftl3.cp file is created, RAID can 
be used as previously discussed to display the flow properties. An anomaly of 
RAIDCONV is that it does not recognize kt=20 type wakes, used where wakes from a 


wing abut against a body [Ref. 4]. To assist in extracting all the data, the kt=20 wakes 


aay 


should be placed at the end of the input file. In general, it is good practice to place all 
wakes at the end of the input file when using A5021. 

The arbitrarily named output file duplicates the data found in the ftl3 file and 
includes moments and forces. A5021 will sum up the moments and forces on each network 
and for all networks so far [App. D]. The moments are computed based on the coordinates 


entered into the input file [Ref. 4and App. A]. 
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IV. RESULTS OF A502i COMPUTATIONS 


A. PARABOLIC ARC AIRFOIL DISCUSSION 


This simple geometry was analyzed primarily to evaluate A502i’s capabilities by a 
comparison to known linear theory. To this end, the geometry discussed in section II-C 
and shown in Figure 2.1 was run by A5021 at a Mach of 0.3 and a Mach of 1.5 at an angle 
of attack of zero. Two of the 49 flow properties that A5021 computes [Ref. 4] for each 
panel are linear Cp and second order Cp, given by 

CPLIN = -2u, (4.1) 

CP2ND = -2u,-[(1-Mx2)u,2+v,.2+W,2] (4.2) 

Where uc, Ve and we are the compressible components of the perturbation velocity. 

Figure 4.1 plots the linear theory, A5021 linear and second order results for the subsonic 
case, while Figure 4.2 represents the supersonic solution. 

Linear theory for parabolic arc airfoils is outlined in Refs. 2 and 3. The equation 


representing the subsonic case is given by: 


—-S*t l-x 
Cp(x) = as (1 -(.5-0) #In |) where 0<x< | (4.3) 


x * chord * ,/ 1 — M~ 





Equation 4.3 includes a Prandtl-Glauert compressibility correction. The equation 


representing the supersonic case 1s given by: 
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Figure 4.2 Cp Comparison of a Parabolic Arc Airfoil at Mach = 1.5. 
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Figure 4.1 shows very good agreement with the linear theory curve. There is a 
small but noticeable difference between A502i’s linear results and second order results, 
with the second order results being more accurate, as expected. The maximum difference 
between the linear theory curve and the second order A502i curve amounts to a value of 
2.5% right at mid-chord. The gap between the two curves from .25 chord to .75 chord 
were the result of thin panel density in that region. 

Figure 4.2 shows excellent agreement with the linear theory curve. The A502! 
values of Cp for the linear and second order analysis are virtually identical. There are small 
deviations from linear theory near the leading and trailing edges, but this is expected due to 
numerical error associated with the discontinuity A5021 would encounter right on the 


leading or trailing edges. 
B. DELTAWING DISCUSSION 


Reference 9 provides data on Mach distribution, using approximated linear theory, 
over a deltawing of the configuration discussed in section II-D. This simple geometry 
provided another test of A502i’s capabilities. Figures 4.3 and 4.5 show the A502i results 
for the subsonic and supersonic case, while Figures 4.4 and 4.6 reflect the results from 
Ref. 9. For both cases, good agreement is found with the linear theory, with A5021’s 
subsonic analysis being physically more accurate than the approximate linear theory, while 
A502i’s supersonic analysis is not as physically accurate. 

A comparison of Figures 4.3 and 4.4 reveals several points of interest. The Mach 
contour representing the free-stream value is given by the dashed line. All lines outside the 
dashed line represent areas where the Mach value is less than free-stream, and inside the 
dashed line is where the Mach value is more than free-stream. The location of where the 
free-stream Mach contour, in Figure 4.3, intersects the centerline agrees very well with 


Figure 4.4. However, Figure 4.4 does not have the contour extending all the way to the 


on 


tip. This is a physical limitation of the approximate theory used in Figure 4.4 and A502i is 
giving a more realistic solution. Figure 4.4 suggests that the peak Mach value occurs at 
approximately two-thirds chord along the centerline. The A502i results show the peak 
Mach contour occurring out midway along the semi-span. Those Mach values are less than 
1% larger than the yellow Mach contour surrounding it, and can be attributed to how the 
panel density increases with movement towards the wingtip. A502i performed very well 
for this subsonic case. 

A companison of Figures 4.5 and 4.6 shows that A502: did not perform as well as 
in the subsonic case. Again, the contour representing the free-stream value of 
Mach =1.414 is given by the dashed line. All lines forward are below free-stream and all 
lines aft are above free-stream. Figure 4.6 shows the intersection of the free-stream Mach 
contour on the centerline occurring at approximately 39% chord, which is in excellent 
agreement with A502i’s result. Figure 4.6 shows the peak Mach value occurring at the 
trailing edge on the centerline. This makes more physical sense than the results 
that A502i yielded. The maximum thickness of the deltawing occurs along the centerline, 
allowing for greater expansion. The discrepancy may be attributable to panel density and 
accumulation of numerical error. A close study of the A502i results reveals some 
discontinuities along the column of panels out at the wing-tip which would have adversely 


affected the solution and caused errors to propagate along the semi-span. 
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Figure 4.3 A502i Mach Contour Plot (M~=.01) 





Figure 4.4 Approximate Linear Theory [Ref. 9] M.=0 
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Lines of Constant Mach (Minf=1.414) 
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Figure 4.5 A502i Mach Contour Plot (M.o=1.414) 
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Figure 4.6 Approximate Linear Theory [Ref. 9] M.=1.414 
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C. GBU-24 FREE-STREAM (NO CANARDS) DISCUSSION 


A free-stream measurement of separation forces on the GBU-24, without canards, 
was conducted tn a wind tunnel for various Mach numbers from .8 to 1.2 [Ref. 7]. Since 
A5021 uses linear potential theory, the model of the GBU-24, without canards, was 
evaluated at both Mach .8 and 1.2, avoiding the transonic regime, to ascertain the accuracy 
of the code with the given geometry. Normal forces and pitching moments for both cases 


are plotted and compared to the wind tunnel data. 


1. Subsonic Case (Mxo=0.8) 


The GBU-24 model, without the canards, was run for angles of attack varying 
from -10 to +10 degrees in two degree increments. Values much higher than that ran into 
wake modelling problems as the wake’s angle relative to the free-stream was getting large 
enough that results would become questionable, and remodelling the wake was too difficult 
for such a complex geometry. The results of the A502: analysis are displayed in Figures 
4.7 and 4.8. For angles of attack between -4 and +4 degrees, A5021 does a good job of 
predicting the separation forces. The pitching moment, which happens to be unstable 
without the canards, is approximately linear over the -4 to +4 degree range and 1s the 
limiting factor to the models accuracy. The normal force is approximately linear over a 
wider range, and A502: does a good job of predicting the normal forces from -6 to +6 


degrees. 


2. Supersonic Case (Mo=1.2) 


The results for the subsonic case demonstrated that the effective range of angle of 
attack that A502i needed to explore was from -6 to +6 degrees. Figures 4.9 and 4.10 plot 
the comparison of wind tunnel data versus A502i results for pitching moment and normal 
force. The results for the supersonic case are slightly better than that of the subsonic case. 


The actual pitching moment of the GBU-24 is approximately linear over a wider angle of 
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Figure 4.8 Comparison of Normal Forces 


36 





attack region, but fluctuations tn the data at -6 and +6 degrees means that the model is still 
only viable from -4 to +4 degrees. The normal force line is nearly linear from -10 to +10 
degrees and extrapolating the A502i results out to 10 degrees would still yield good 


predictions. 
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“Figure 4.10 Comparison of Normal Forces ; 
D. GBU-24 FREE-STREAM (WITH CANARDS) DISCUSSION 

As in the case with no canards, free-stream measurements of the separation forces 
on GBU-24 were taken from Mach .8 to 1.2 in a wind tunnel [Ref. 7]. Again, due to the 
limitations of linear theory inherent in the code, an analysis was done for Mach numbers of 
8 and 1.2 to minimize transonic effects. Even with the more complex geometry, A5021 
does an accurate job of predicting the separation forces over the range of angles of attack 


that are approximately linear. 


1. Subsonic Case (Mo=0.8) 


The GBU-24 model, with canards, was run in two degree increments of angle of 
attack from -10 to +10. The wake modelling limitation, as well as reviewing the data from 


the wind tunnel measurements [Ref. 7] showed the non-linearity of the separation forces at 


38 





the higher values of angle of attack, precluded any attempts to predict forces beyond the 


aforementioned angle of attack interval. The results of the A502i analysis are displayed in 


Figures 4.11 and 4.12. The addition of the canards makes the pitching moment stable, but 


linear over a smaller region than without the canards. A5021 gave accurate results from -2 


to +2 degrees angle of attack when predicting pitching moment. The prediction of normal 


forces fared better, showing accurate results from -3 to +3 degrees angle of attack. 
GBU-24 Free-Stream (M=0.8) 
Canards On 
Pitching Moment 





—@— Wind Tunnel Dat: 
—#— A502 Results 


Alpha (Deg) 


Figure 4.11 Comparison of Pitching Moments 
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Figure 4.12 Comparison of Normal Forces 


2. Supersonic Case (Mo=1.2) 
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As in the case with no canards, the region of accuracy, with the model used, was 


assumed to be less than + or - 10 degrees angle of attack. Cases were run from -6 to +6 


degrees angle of attack in two degree increments. A comparison of A502i results with 


wind tunnel data is shown in Figures 4.13 and 4.14. For both the pitching moment and the 


normal force, A502i does a much better job of prediction then when subsonic. The wind 


tunnel data is nearly linear in both pitch moment and normal force from -8 to +8 degrees 


angle of attack. Extrapolating the A502i data out to + or - 8 degrees angle of attack, shows 


excellent agreement with the wind tunnel data. 
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Figure 4.14 Comparison of Normal Forces 
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E. F-14 DISCUSSION 


The GBU-24, without canards was located at station 3 of the F-14 via the NAVSEP 
code. In the process of running a solution, the combined geometries were found to have a 
total of 150 networks (A502i1's maximum). The combined geometry had about 4,000 
panels, far short of the 20,000 panel maximum, so there was room for more detailed 
modelling, but there was not a chance to insert the GBU-24, with canards, into the F-14’s 
flow field and analyze the forces on the bomb. The geometry with the canardless bomb ran 
to a solution that appeared to be valid, so there 1s a high degree of confidence that if the 
number of networks could be reduced to allow the GBU-24, with canards, to be inserted 
into the F-14 flow field, the code would yield accurate predictions at small angles of attack 
on the forward stations. To reduce the total number of networks by combining existing 
networks would have required a large time investment and the use of MACGS, which the 
department currently does not possess, the two reasons why it was not done. Figure 4.15 
shows a Mach distribution of the solution of the canardless bomb and F-14 at Mach = 0.8 


and O degrees angle of attack 
F. POST-PROCESSING DISCUSSION 


The Mach values for the subsonic case of the GBU-24, with canards, at 4 degrees 
angle of attack, were extracted from the ft13 file. These values, used in conjunction with 
RAID are shown in Figures 4.16 and 4.17. The color distribution over the nose in Figure 
4.16 indicates that the bomb is at an angle of attack, and scanning the rest of the model 
showed no discontinuous solutions, which is generally represented in A5021 by a Mach 
value of O or 1,000. The visual representation is a quick way of telling if A502i ran an 
accurate solution. The only other way is to individually check the Mach or Cp values of 
each panel in the ft13 file or the arbitrarily named output file. The other point of interest in 


Figure 4.16 is the lack of panel density along the mid-section of the bomb. The goal, in the 
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case of stores separation prediction, ts to have as simple a model as possible that still gives 
accurate predictions. The fewer the number of panels, the shorter the run time. The fact 
that A502 is a higher order panel method allows the luxury of using fewer panels. Figure 
4.17 highlights the approach used to take into account separation effects as discussed on 


page 16. 


43 





) rts 4. 
ae a 


Figure 4.15 Mach Distribution over GBU-24 and F-14 (M.=0.8, 
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4.17 Mach Distribution over GBU-24 (M.=0.8, alpha=4°) 








V. SUMMARY AND CONCLUSIONS 


The main goal of this analysis is to determine the accuracy of A502i on both simple 
geometries and complex geometries. To accomplish this, A502i is compared with results 
available from linear theory and wind tunnel experiments. This allows conclusions to be 
drawn on the capabilities as well as the limitations of A5021. 

In general, A502i can eeeuraicly predict flow properties, forces and moments on 
simple and complex geometries at low angles of attack. The predictions are valid over a 
wide Mach range, from 0 up to and including 0.8 and from 1.2 and above. The supersonic 
solutions are available due to A5021’s higher order capabilities. 

The limitations of A5021 are consistent with most panel methods. A502i cannot 
predict flow dominated by viscous, separated or transonic effects. It cannot predict flow 
with different total pressures, such as flow properties inside a jet plume or a propeller 
slipstream swirl. The biggest shortcoming of A5021 is its inability to handle unsteady cases 
and automatically determine wake shapes. 

Experience or knowledge of the flow properties around the geometry being tested is 
important in building an accurate model. An accurate model may not be physically 
accurate. Flight test results revealed a yawing moment on the GBU-24 that was not 
discovered in wind tunnel experiments when the bomb was carried on an aft station. The 
yawing moment may be caused by the fact that the canards are not fixed, but spring- 
damped. A502i predictions would not be accurate without inputting a moment to simulate 
the deflection of the canards, since the canards are fixed by the geometry. A5021, while a 
powerful tool in terms of cost savings and time, cannot completely substitute for wind 
tunnel experiments and flight tests, as constructing a complex geometry to exact physical 


specifications will probably not yield accurate predictions. 
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APPENDIX A. GBU-24 OUPUT FILE (INPUT FILE PORTION) 
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dblt edge starts at ai # 31 ends at ai # 29 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
162 eZ 4 1 2 3 4 
25.4 12 4 7 6 5 4 
3251 18 5 mu-match 7 8 9 10 
Oabutment # 47 
dblt edge starts at ai #* 32 ends at ai # 26 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
ey on ie 4 2 3 4 -5 -6 a 8 Sy 10. oil 
18.4 12 4 11 -10 9 8 -7 -6 5 4 3 z aL 
Zi i 4 mu-match Lie 1) 9 8 -7 -6 S 4 3 2 a 
Oabutment # 48 
dblt edge starts at ai # 32 ends at ai # 31 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
Ut ae 12 4 iL 2 3 4 
2354 12 4 10 9 8 7 
a2 ak 18 5 mu-match 4 5 6 7 
Oabutment # 49 doublet strength matched to zero along this abutment ae Warning ** 
dblt edge starts at ai # 26 ends at ai # 33 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
18 4% 12 4 mu-match 1 2 3 4 s 6 7 
OQabutment # 50 doublet strength matched to zero along this abutment = warning ** 
dblt edge Starts at al # 33 ends at ai # 34 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
18a2 12 4 mu-match 1 2 3 4 5 6 7 8 9° 10 Ta 
Qabutment # Sl 
adblt edge starts at ai # 32 ends at ai # 34 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
18.3 12 4 vor-mtch 1 2 3 4 5 6 7 
2¢ 74 18 5 mu-match vor-mtch a 6 5 4 3 2 2 
Oabutment # 52 doublet strength matched to zero along this abutment "““sewarning ** 
dblt edge starts at ai # 28 ends at ai # 35 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
13. y 12 4 mu-match 1 2 3 4 5 6 2 
Oabutment # 53 doublet strength matched to zero along this abutment —=—eevarning  ** 
dblt edge starts at ai # 35 ends at ai # 36 
nw.edge nw/id type type matching kutta-fl corresponding edge points ({ minus (-) indicates point moved by $eat 
Toe i 4 mu-match ah 2 3 4 5 6 q 8 9 10 22 
Oabutment # 54 
dblt edge starts at ai # 29 ends at ai # 36 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
Lg. 3 a2 4 vor-mtch i 2 3 4 5 6 7 
2620 18 5 mu-match vor-mtch a 6 5 4 3 2 i 
Oabutment # 55 doublet strength matched to zero along this abutment _ eewracning. ** 
dblt edge Starts at ai # 27 ends at ai # 37 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
20 al 12 4 mu-match I: 2 3 4 5 6 fi 
Oabutment # 56 doublet strength matched to zero along this abutment *** warning ** 
* 
dblt edge starts at ai # 37 ends at ai # 38 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
20:22 Zz 4 mu-match i 2 3 4 5 6 7 8 9-10 11 
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dblt edge starts at ai # 45 ends at ai # 30 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
25.4 18 2 ao 2 
S22 18 2 2 ¥ 
S25 18 2 1 2 
Oabutment # 69 wake side edge left unabutted [Fe warning ** 
dblt edge starts at ai # 36 ends at ai # 46 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
26.2 18 Z sf iz 
Oabutment # 70 wake trailing edge unabutted. wake filaments will be added *** gentle reminder ** 
dblt edge Starts at ai # 46 ends at ai # 47 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
26.3 18 Z il Pe 3 4 5 6 7 
Oabutment # 71 wake side edge left unabutted *** warning ** 
dblt edge starts at ai # 47 ends at ai # 29 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
26.04 18 z 1 2 
Oabutment # 72 wake side edge left unabutted 7 warning ** 
dblt edge starts at ai # 38 ends at ai # 48 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
2 hie 18 2 1 2 
Oabutment # 73 wake trailing edge unabutted. wake filaments will be added *** gentle reminder ** 
dblt edge starts at ai # 48 ends at ai # 49 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
27.3 18 2 a 2 3 4 5 6 7 
Oabutment # 74 wake side edge left unabutted fe varning ** 
dblt edge starts at ai # 49 ends at ai # 31 
nw.edge nw/id type type matching kutta-fl corresponding edge points ({ minus (-) indicates point moved by Seat 
27.4 18 2 1 2 
Oabutment # 75 wake side edge left unabutted warning *~* 
dblt edge starts at ai # 14 ends at ai # 50 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
28.2 18 2 i 2 3 4 5 6 a 8 9eerO. 11 25 13 14 #15 "T6 
Oabutment # 76 wake trailing edge unabutted. wake filaments will be added *** gentle reminder ** 
dblt edge starts at ai # 50 ends at ai # 51 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
2607s 18 2 1 2 3 4 5 6 7 8 9 
Oabutment # 77 wake side edge left unabutted =z warning —** 
dblt edge starts at ai # 51 ends at ai # 13 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by Seat 
28.4 18 2 1 Zz 2 4 5 6 7 8 Ss 10 Jy 12 35 14 15 is 
Oabutment # 78 wake side edge left unabutted [= Swarning ** 
dblt edge starts at ai # 6 ends at ai # 52 
nw.edge nw/id type type matching kutta-fl corresponding edge points ( minus (-) indicates point moved by $eat 


Zoe 18 Z 1 2 S 4 5 6 a 8 Seno “bi 12 is 14 “15 75t6 
Oabutment # 79 wake trailing edge unabutted. wake filaments will be added *** gentle reminder ** 
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APPENDIX C. GBU-24 OUTPUT FILE (UNIT NORMALS) 


6°O & 


oO - 


ej Seszy janqe opzt usf uzyt oz ds‘ teued mu 


91, 88eg 


v 


i 


ee 


ce 


4 


cc 


T 
i 


Le 


TE 


Aeidso wo Jnojeyew Aq peqwuid 


LT 
ot 
vT 
cL 


OT 


c 
+ 


ot 
eT 
EL 
GL 
Lt 


OT 


© 
ie 


Of 


Of 


T 


T 


6c 


6c 


ElcLe. Oo 
600L9L'°0 
6POL9L‘O 
LE69924°0 
LCCI9L 6 
CSSLIOL‘O 
CSSL9L'O 
660L9L°0 
EGELIL 6 
EOELSL 0 
ShcEese, 0 
600L9L°0 
600L9L°0 
EC6SIE 0 


G67994 70 


CGELUS 0 
88STh9°O 
88STP9°O 
VELTIS 0 
PELTYS 29 
L86009'0 
L86009°0 
OSOTP9'O 
SOcLeS 0 
BETTH9'O 
T6ETHI'O 
88STP9'O 
88STP9°O 
VELUpSo 6 


8B8Th9°O 


ouz 


8c 


Bz 


6¢T000°0 
000000°0 
000000°0 
6¢T000°0 
6¢T000°0 
000000°0 
000000°0 
6¢T000'0 
000000°0 
000000°0 
6¢T000°0 
000000°0 
000000°0 
6¢T000°0 


000000°0 


Le 


Le 


860606 °C- 
680606 °C- 
OSb6R6°Z- 
£89606 °C- 
VELEC6 Ze 
TODCOb Z- 
STCUOr <= 
STSZOP Z- 
Cosccur cz 
CISCOP <c- 
O8ScOP’°c- 
ScoCcOP c- 
Sc9CcOP <c- 
PRSCOr c- 


LLGUT? ¢c- 


92 


a6 


OL8SB88°E GS96TC'9E ST ST 
000 
tosses’ ec B8TPEB8L' SE PT PT 
000 
9SLS88°E CPLESE'TP €T cal 
000 
9ELSBB°E LCTOEG ED CL cL 
Sc9 
S608S8°E¢ 96°S80°SP EE ie 
000 
6CELEC & STPETO’ OS OT OT 
S29 
COPOET'E LECeso Se 6 6 
S248 
O00cET EC SLUTLEL Se 8 8 
S2T 
oo00cTET' EC O6S6RP'OL & c 
SLt 
OOOCET'E COLLEGE FE 3 9 
S29 
OO0O0TET EC beSbos'SE S S 
SL8 
O0O0CET EC TTST8SP' BE P v 
SZT 
OOOCET'E COPSST TP € te 
SLE 
oodcwE? EC 8T8PEs EP c c 
000 
ELOOPT'E S6€060'SP ac ic 
n 
22 eAteus/of oc 
T >: YAOMAEU TOJ squtod Torzquos Oo 
qeoehqtt.®.0 
~£0000'°0 Tue/Teojqe 
Sct660'°C upt3qe/ 3 
ce 
ee 
O€ 
62 
82 
bc 
Se 
$2 
PZ 
Sé bz ce 


stsAjTeue sue{td zajjez, OJ syIOMAeU Jo Areuuns 


ae 


§NO/JAAT MU 


— ee ee ee ee —_——_—— 


PT-yx2OmJeuU MU 


bz a 
uT/JTATMuU 





LSE SBEh 2 JEW 





85 


86 





vr 
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